Observation of nonadditive mixed state phases with 
polarized neutrons 
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In a neutron polarimetry experiment the mixed state relative phases between spin eigenstates are 
determined from the maxima and minima of measured intensity oscillations. We consider evolutions 
leading to purely geometric, purely dynamical and combined phases. It is experimentally demon- 
strated that the sum of the individually determined geometric and dynamical phases is not equal 
to the associated total phase which is obtained from a single measurement, unless the system is in 
a pure state. 
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Evolving quantum systems acquire two kinds of phase 
factors: (i) the dynamical phase which depends on the 
dynamical properties of the system - like energy or time 
- during a particular evolution, and (ii) the geometric 
phase which only depends on the path the system takes 
in state space on its way from the initial to the final 
state. Since its discovery by Pancharatnam [l| and Berry 
the concept was widely expanded and has undergone 
several generalizations. Nonadiabatic [3| and noncyclic 
[3 evolutions as well as the off-diagonal case, where ini- 
tial and final state are mutually orthogonal [3| , have been 
considered. Ever since, a great variety of experimental 
demonstrations has been accomplished [1, 0| , also in neu- 
tron optics [l,[l,[l0l. [Tl|. Due to its potential robustness 
against noise |12l | the geometric phase is an excellent can- 
didate to be utilized for logic gate operations in quantum 
information science p^ . Thus, a rigorous investigation 
of all its properties is of great importance. 

In addition to an approach by Uhlmann a new 
concept of phase for mixed input states based on inter- 
ferometry was developed by Sjoqvist et al. 15]. Here, 
each eigenvector of the initial density matrix indepen- 
dently acquires a geometric phase. The total mixed state 
phase is a weighted average of the individual phase fac- 
tors. This concept is of great significance for such exper- 
imental situations or technical applications where pure 
state theories may imply strong idealizations. Theoret- 
ical predictions have been tested by Du et al. [lB| and 
Ericsson et al. [13] using NMR and single-photon in- 
terferometry, respectively. In this letter, we report on 
measurements of nonadiabatic and noncyclic geometric, 
dynamical and combined phases. These depend on noise 
strength in state preparation, defining the degree of po- 
larization, the purity, of the neutron input state. In 
particular, our experiment demonstrates that the geo- 
metric and dynamical mixed state phases $g and <I>d, 
resulting from separate measurements, are not additive 
p^ . because the phase resulting from a single, cumula- 
tive, measurement differs from $g -I- $d- This nonaddi- 
tivity might be of practical importance for applications, 
e.g. the design of phase gates for quantum computation. 

A neutron beam propagating in y-direction interact- 




FIG. 1: Sketch of neutron polarimetry setup for phase mea- 
surement with overall guide field Bz, polarizer P, three DC- 
coils to implement unitary operations Ui, U|, U^, analyzer 
A and detector D. Greek letters denote spin rotation angles. 
Shifting the second coil induces an additional dynamical phase 
rj resulting in intensity oscillations. The desired phase (f> is de- 
termined from their minima and maxima. 

ing with static magnetic fields B(y) is described by the 
Hamiltonian H= —h^/2mW^ — ii<jB{y)] m and /i are 
the mass and the magnetic moment of the neutron, re- 
spectively. Zeeman splitting within B(j/) leads to so- 
lutions of the Schrodinger equation cos(6'/2)]fc+)]+) -I- 
e'" sin(0/2)]fc_)]— ), where \k±) are the momentum and 
]±) the spin eigenstates within the field B(?/). 9 and a 
denote the polar and azimuthal angles determining the 
direction of the spin with respect to B(y). k± ~ fco T 
where fco is the momentum of the free particle and Ak — 
m^\B{y)\/h^k is the field-induced momentum shift. Ak 
can be detected from spinor precession. Omitting the 
coupling of momentum and spin, we focus on the evolu- 
tion of superposed spin eigenstates resulting in Larmor 
precession of the polarization vector r = {ip\a\ip), where 
(? — {cxjO-y^az) is the Pauli vector operator. The uni- 
tary, unimodular operator 

U(C',<5',C') = e^''co8^'\+){+\-e~^<'sme\+)H 
+ e< sin e' 1 -) (+ 1 + e-'^' cos C' ] - > (- 1 
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FIG. 2: Evolution of the |+) eigenstate induced by U^, asso- 
ciated to: a) Purely (noncyclic) geometric phase (2^ — it/ 2). 
b) Combinations of dynamical and geometric phase on the 
Bloch sphere (0 < 2^ < 7r/2). 



describes the evolution of the system vifithin static mag- 
netic fields. The set of SU(2) parameters {£,',S',C) is 
related to the so-called Cayley-Klein parameters a, b via 
a = e"' cos^' and b = — e""" sin^' (see e.g. [l9|). 

Consider the experimental setup shown in Fig.[TJ A 
monochromatic neutron beam passes the polarizer P 
preparing it in the 'up' state |-|-) with respect to a mag- 
netic guide field in z-direction (B^). Next, the beam 
approaches a DC coil with its field B^; pointing to the 
x-direction. Bj; is chosen such that it carries out the 
transformation Ui = U(7r/4, 0, — 7r/2), corresponding to a 
+Tr/2 rotation around the +a;-axis. After Ui the result- 
ing state of the system is a coherent superposition of the 
two orthogonal spin eigenstates: \ipo) = l/^/2{\+)—i\—)). 
A subsequent coil, represented by U(^, 0, —tt/2), is set to 
cause a spin rotation around the -|-a;-axis by an angle 2^. 
This second coil and the following propagation distance 
within B^ - corresponding to a rotation angle 26 around 



the -|-z-axis - define an evolution = U(^, d, Q. Under- 
going the transformation the two spin eigenstates |±) 
acquire opposite total phase ±0 — ± arg(±|U0|±) — ±S. 
A third coil (u|) reverses the action of the first one and 
would therefore transform a state lipo) back to |-|-). The 
state of the system entering the third coil equals \tpo) 
only if U0 1. can be extracted by applying an extra 
dynamical phase shift i^ry to |±). It is implemented by 
adjusting both inter-coil distances from first to second 
and second to third coil to spin rotation angle equiva- 
lents of 27rn + T] and 27rn' + 26 — 1] (n, n' are integer). 
By scanning the position of the second coil these rota- 
tion angles - referred to as U,, = U(0, 77/2) and U]^ - are 
varied to yield intensity oscillations from which <j) is cal- 
culated. Note that, because of ^' = 0, the parameter 
C' is undetermined and, therefore, omitted in U^. Af- 
ter projection on the state j-f) by the analyzer A, the 
phase (j)~5 and its visibility v — \ cos^| can be computed 
as functions of the maxima and minima of the intensity, 
Imax and Imin, mcasurcd by the detector D [iO]. More 
general, a neutron beam with incident purity Tq = |r^| 
along the -|-z-axis (fp — (0,0, Tq)) is described by the 
density operator pin(7'o) = I / 2{\ + r'^a z) ■ Calculating 
Pout = Ulu|,U0U„Uipi„UluJ,U^U^Ui, we find the in- 
tensity 

1 - r' 

1^ (X -Y^+r'a{cos^^cos'^5 + siTi^^cos^{C~v)X^) 

after the analyzer A. Considering the maxima and min- 
ima IJ!,„^, Ij^i^ of n-induced oscillations of F one obtains 

iiidx ^ niiii ' 

the mixed state phase [2l| 



^(,' ) _ arecos . [Ifnin/Ig- 1/2(1 -r(,)]A-- 



with a normalization factor I(J = 2Ig/(l + Tq). Iq is the 
intensity measured at = 1. 

The noncyclic geometric phase is given by cjjg — —Q/2, 
where ft is the solid angle enclosed by a geodesic path 
and its shortest geodesic closure on the Bloch sphere [11 : 
(j>g and the total phase 4> are related to the path by the 
polar and azimuthal angles 2^ and 25, so that the pure 
state geometric phase becomes 

0g = <^-(/)d = <5[l-cos(2O], (3) 

while its dynamical counterpart is 

0d = (5cos(2O. (4) 

By proper choice of 2^ and 2(5, can be set to generate 
purely geometric, purely dynamical, or arbitrary combi- 
nations of both phases as shown in Fig. [51 



The theoretical prediction for the mixed state phase is 

MM 

$(rQ) = arctan (rQ tan (5) (5) 

To access Eq. ([S|) experimentally rg needs to be varied. 
In addition to the DC current, which effects the trans- 
formation Ui, random noise is applied to the first coil, 
thereby changing B^, in time. Neutrons, which are part 
of the ensemble pinifo), arrive at different times at the 
coil and experience different magnetic field strengths. 
This is equivalent to applying different unitary operators 
U(7r/4 + AC'(i),0,-7r/2) = Ui(A,f'(i)). For the whole 
ensemble we have to take the time integral 

p = J Ui(Ae'(i))|+)(+|uI(A^'(i))di. 
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FIG. 3: Above: typical measured intensities for five values 
of purity ro versus coil positions. Below: resulting mixed 
state phases <l?g and versus purity tq. a) Purely geometric 
case: 2^ = 7r/2. Data indicated by circles is calculated from 
oscillations above, b) Purely dynamical case: 2^ — 0. Only 
one intensity value - marked by rectangles in the above plot 
- is needed for each $jj (see text). 



Although each separate transformation is unitary, due to 
the randomness of the signal we end up with a nonuni- 
tary evolution that yields a mixed state Note that 
in this method the purity t'q of the input state is not af- 
fected before Ui creates spin superpositions distributed 
around ItJjq) within the y, z-plane. We are left with 
ro = (0, — ro, 0) where rg < 1, as has been confirmed by a 

3D spin analysis of the state UipinUi. 

The experiment was carried out at the research reac- 
tor facility of the Vienna University of Technology. A 
neutron beam ~ incident from a pyrolytic graphite crys- 
tal - with a mean wavelength A ~ 1.98 A and spectral 
width AA/A ~ 0.015, was polarized up to rg « 99% by 
reflection from a bent Co-Ti supermirror array. The an- 
alyzing supermirror was slightly de-adjusted to higher in- 
cident angles to suppress second order harmonics in the 
incident beam. The final maximum intensity was about 
150counts/s at a beam cross-section of roughly 1 cm^. A 
•^He gas detector with high efficiency for neutrons of this 
energy range was used. Spin rotations around the +x- 
axis were implemented by magnetic fields B^; from DC 
coils made of anodized aluminium wire (0.6 mm in diam- 
eter) wound on frames with rectangular profile (7x7x2 
cm'^). Coil windings in ^-direction provided for compen- 
sation of the guide field at coil positions (see Fig.[T|) 
for Bj; to be parallel to the x-direction. was realized 
by two rectangular coils of 150 cm length in Helmholtz 
geometry. Maximum coil currents of 2 A and guide field 
strengths of 1 mT were sufficient to achieve required spin 



rotations and prevent unwanted depolarization. The 
noise from a standard signal generator consisted of ran- 
dom DC offsets varying at a rate of 20 kHz. In order 
to find the coil current values for required spin rotation 
angles, each coil current was scanned separately. A mea- 
sured intensity minimum after stepwise increase indicates 
that a TT rotation has been achieved. Then, with both 
coils mounted on translation stages and set to a 7r/2 ro- 
tation around the a;-axis, inter-coil distances were varied 
to search for intensity minima. The distance between two 
minima (6-7 cm in our case) corresponds to a 27r rotation 
around B^ within the x, y plane. 

For purely geometric phases the parameter sets (^ — 
7r/4, (5, C = (5 - 7r/2) with 5 = 0g = 7r/8, 27r/8, Stt/S were 
chosen. For each set the intensity oscillations V (see 
Fig. [3^, upper graph; data shown for J = 0g = 7r/8) were 
measured, scanning the position of the second coil for 
five values of rg . For the purely dynamical phase - since 
2^ = and Eq. ^ reduces to V ^ {l-r'^) /2+r'Q cos^ 5 - 
only one intensity value is needed. With the second coil 
turned off, the distance between the first and the third 
coil was chosen such that 5 — (fi^i — 7r/8, 27r/8, 37r/8 for 
five values of rg (Fig.[3jD, upper graph). By Eq. ^ the 
geometric and dynamical mixed state phases 5'g(rg) and 
^d{ro) (Fig. [21 lower graphs) were calculated from inten- 
sity values extracted from least square fits for V (solid 
lines in Fig. [31 upper graphs). Vertical error bars in the 
lower graphs contain two about equal contributions: fit- 
ting errors and an estimated uncertainty of 0.5 mm for 
the reproduction of coil positions. The small horizontal 
error bars stem from the statistical uncertainty of the pu- 
rity measurements. The solid lines are theoretical curves 
according to Eq. ([5]), using the measured value for <i> with- 
out noise as phase reference. The experimental data 
reproduce well the rg-dependence predicted by Eq. ([5]). 

We want to emphasize that our investigation focuses 
on a special property of the mixed state phase: its 
nonadditivity. Since the Sjoqvist phase is defined as 
a weighted average of phase factors rather then one of 
phases (see [3, for a more elaborate discussion) it 
is true only for pure states that phases of separate mea- 
surements can be added up to a total phase. Suppose 
we carry out two measurements on a pure state sys- 
tem: the state is subjected to a unitary transforma- 
tion Ug in the first and to Ud in the second measure- 
ment inducing the phases 4>g and 4>d, respectively. Ap- 
plying Eqs. ([3]) and (|4]), we can also choose a combina- 
tion of angles 2^ and 25 leading to a transformation Utot 
so that we measure the total pure state phase (j)^ -\- 4>d 
(note that the three evolution paths induced by Ug, Ud 
and Utot differ from each other). However, the result 
of the latter experiment for the system in a mixed state 
is $tot(?'o) — arctan [rg tan(0g -I- (/id)]. The total phase 
is then not given by $g(?'g) -I- ^difo), with $g(ro) = 
arctan (rg tan0g) and ^diro) = arctan (rg tan 0d). In 
our experiment we have chosen two examples of Utot, 
i.e. two sets of values for B^, in the second coil (2f '^^^ = 
60°, 2^(2) ^ 480-) ^Yie distance within B^ {26^^^ = 
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FIG. 4: Filled markers: Measured total mixed state phase 
$tot versus purity ro for two examples of Utot associated to 
the total pure state phases 0g + "^J^ (see text). Open 

markers: $g + $j ' as calculated from measured data in 
Fig. |3] The solid and dotted lines are theory curves assum- 
ing either nonadditivity or additivity, respectively. The data 
clearly demonstrate the nonadditivity of the Sjoqvist mixed 
state phase. 

90°, 25(2) = 135°). According to Eqs. © and dl the 
total pure state phases ^^^^ + (\)^^ and 0g^^ + <i)f^ with 

(\)^^^ — — 7r/8 and ~ 27r/8, 7r/8 are obtained 

from intensity oscillations. In Fig.|4] the resulting mixed 
state phases <^\^^ and the sum are plotted. 

Note that this nonadditivity of mixed state phases is 



not due to the nonlinearity of the geometric phase, that 
occurs - for instance - when the system evolves close to 
the orthogonal state of the reference state [l^ • 

Recently there has been a report on NMR experi- 
ments [25j investigating Uhlmann's mixed state geomet- 
ric phase. It is a property of a composite system under- 
going a certain non-local evolution of system and ancilla 
[26j . Diverse phase definitions, depending on this evolu- 
tion, are possible. The phase investigated in the present 
paper is a special case in which the ancilla does not 
necessarily evolve. While the preconditions for inherent 
fault tolerance remain intact for the Sjoqvist phase, the 
question whether other phases offer advantages in terms 
of robustness remains an exciting issue of discussion. 

To summarize, we have measured spin-1/2 mixed state 
phases with polarized neutrons. Their dependence on 
the purity of the input state was observed. Our results 
show that these phases are not additive. The sum of 
phases measured in separate experiments is not equal 
to the result of one single measurement, in contrast to 
what could be expected from straightforward extrapo- 
lation of pure state behavior. This interesting prop- 
erty of the Sjoqvist mixed state phase might be of high 
relevance for possible applications of geometric phases. 

We thank E. Balcar for critical reading of the 
manuscript. This work was financed by the Japan Sci- 
ence and Technology Agency (JST) and the Austrian Sci- 
ence Fund (FWF, Project. Nr. P 17803-N02). 



[1] Pancharatnam, Proc. Indian Acad. Sci. A 44 (1956). 
[2] M.V. Berry, Proc. R. Soc. Lond. A 392, 45 (1984). 
[3] Y. Aharonov and J.S. Anandan, Phys. Rev. Lett. 58, 
1593 (1987). 

[4] J. Samuel and R. Bhandari, Phys. Rev. Lett. 60, 2339 
(1988). 

[5] N. Manini and F. Pistolesi, Phys. Rev. Lett. 85, 3076 
(2000). 

[6] A. Tomita and R.Y. Chiao, Phys. Rev. Lett. 57, 937 

(1986) . 

[7] D. Suter, K.T. Mueller and A. Pines, Phys. Rev. Lett. 

60, 1218 (1988). 
[8] T. Bitter and D. Dubbers, Phys. Rev. Lett. 59, 251 

(1987) . 

[9] B.E. AUman et al, Phys. Rev. A 56, 4420 (1997). 
[10] Y. Hasegawa et al, Phys. Rev. Lett. 87, 070401 (2001). 
[11] S. Filipp, Y. Hasegawa, R. Loidl and H. Ranch, Phys. 

Rev. A 72, 021602(R) (2005). 
[12] P.J. Leek et al, Science 318, 1889 (2007). 
[13] M.A. Nielsen and I.L. Chuang, Quantum Computation 

and Quantum Information, Cambridge University Press, 

Cambridge, England, 2000. 



[14] A. Uhlmann, Lett. Math. Phys. 21, 229 (1991). 

[15] E. Sjoqvist et al., Phys. Rev. Lett. 85, 2845 (2000). 

[16] J. Du et al, Phys. Rev. Lett. 91, 100403 (2003). 

[17] M. Ericsson et al, Phys. Rev. Lett. 94, 050401 (2005). 

[18] K. Singh et al, Phys. Rev. A 67, 032106 (2003). 

[19] J.J. Sakurai, Modern Quantum Mechanics (Addison- 

Wesley, New York, 1994) 
[20] A.G. Wagh and V.C.Rakhecha, Phys. Lett. A 197, 112 

(1995). 

[21] P. Larsson and E. Sjoqvist, Phys. Lett. A 315, 12 (2003); 
J. Klepp et al, Phys. Lett. A 342, 48 (2005); S. Sponar 
et al.. Acta Phys. Hung. A 26, 165 (2006). 

[22] R.A. Bertlmann, K. Durstberger and Y. Hasegawa, Phys. 
Rev. A 73, 022111 (2006) 

[23] Li-Bin Fu and Jing-Ling Chen, J. Phys. A: Math. Gen. 
37, 3699 (2004). E. Sjoqvist, J. Phys. A: Math. Gen. 37, 
7393 (2004). Li-Bin Fu and Jing-Ling Chen, J. Phys. A: 
Math. Gen. 37, 7395 (2004). 

[24] R. Bhandari , Phys. R ep. 281, 1 (1997). 

[25] J. Du et aZ.. larXiv:0f l0.5804>^l [quant-ph]. 

[26] M. Ericsson et al, Phys. Rev. Lett. 91, 090405 (2003). 



